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Method for determining the dose of at least one active agent based on a genetic analysis. The method 
comprises analyzing specific genes for their nucleotide sequence or their expression levels of gene- 
specific proteins and/or RNA molecules. The gene-specific data is assigned one or more relevant 
physiological functions of the human or animal body, in particular those which have an influence on the 
metabolism, absorption, excretion or distribution of the active agent in the body. The gene-specific data 
and the assigned physiological functions are then integrated into a physiology-based pharmacokinetic 
model (PBPK model). The PBPK model integrates pharmacokinetic data relating to one or more active 
agents. The PBPK model may also receive and evaluate patient-specific data directly inputted and 
combined with data from a knowledge database comprising known values of pharmacokinetic 
parameters. The PBPK's integration of these data provide a calculation of the individual dose of the active 
agent. 
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(54) Bezeichnung: Verfahren zur Bestimmung einer Wirkstoffdosierung 

(57) Zusammenfassung: Es wird ein Verfahren zur Bestim- 
mung der Dosierung von mindestens einem Wirkstoff auf 
Basis einer genetischen Analyse beschrieben. Das Verfah- 
ren umfasst die folgenden Schritte: 

a) Analyse (101) von spezifischen Gensequenzen, mittels 
eines gensequenz-spezifischen Analysengerates, insbe- 
sondere eines sequenz-spezifischen Sensors, oder Be- 
stimmung der Expression von Proteinen entweder iiber 
RNA Transkription mittels quantitativer RNA-spezifischer 
Nachweismethoden oder direkte Messeung der Proteinex- 
pression durch ein Protein-Analysegerat, 

b) Zuordnung der Gensequenzen zu physiologischen 
Funktionen des menschlichen oder tierischen Korpers, ins- 
besondere solchen physiologischen Funktionen, die Ein- 
fluss auf den Abbau, auf die Aufnahme, auf die Abgabe 
oder auf das Verteilungsverhalten des Wirkstoffes im Kor- 
per haben, 

c) Weitergabe der genetischen und der Zuordnungsdaten 
an ein physiologie-basiertes pharmakokinetisches Modell 
(PBPK-Modell)(108), 

d) Eingabe von Wirkstoff-spezifischen Daten in das PB- 
PK-Modell (108), 

e) Eingabe von charakteristischen Patientendaten, gege- 
benenfalls aus direkten Messungen am Korper, 

f) Berechnung von fur das PBPK-Modell notwendigen phy- 
siologischen Einflussparametern aus den Patientendaten 
unter Verwendung von in der Wissensdatenbank enthalte- 
nen Informationen und Weitergabe der Parameter an das 
PBPK-Modell (108), 

g) Berechnung der individuellen Dosis aus den Daten nach 
den Schritten c), d) und f) unter Verwendung des PB- 
PK-Modells(108). 
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Beschreibung 

[0001] Die Erfindung betrifft ein Verfahren zur Bestimmiing der individuellen Dosierung von Arzneimitteln, fur 
die bekannt ist, dass ihre Wirkung durch eine von individuellen Faktoren der Patienten abhangige Pharmako- 
kinetik und/oder Pharmakodynamik beeinflusst wird. Das Verfahren kann je nach Ausfuhrung sowohl als Point 
of Care Losung direkt in der Klinik oder Arztpraxis oder als spezielle Methode im Labormedizinischen Bereich 
zum Einsatz kommen. 

Stand der Technik 

[0002] Die therapeutische Wirkung von Medikamenten wird sowohl von der intrinsischen, biochemischen Wir- 
kung des Wirkstoffs direkt am biologischen Zielmolekul, als auch von der Konzentration am Wirkort bestimmt. 
Die Konzentration am Wirkort wird wiederum von verschiedenen Faktoren wie absorbierter Anteil bei oraler 
Verabreichung, Verteilung im Korper und Geschwindigkeit des metabolischen Abbaus und der Ausscheidung 
beeinflusst. Diese Prozesse hangen stark von physiologischen und anatomischen Eigenschaften des Korpers 
des behandelten Patienten ab. Im einzelnen sind dabei zu nennen: 

- Volumina und Fettgehalt der einzelnen Organe oder als zusammengefasster Parameter das Korperge- 
wicht und der Korperfettgehalt. 

- Blutflussraten in den einzelnen Organen. 

- Funktion des Magen-Darm-Traktes 

- Funktion der Ausscheidungsorgane wie Niere wie Galientrakt. 

- Expression und Funktion von Abbauenzymen, insbesondere in Leber und Darm 

- Expression und Funktion von Proteinen fur den aktiven Transport von Molekulen durch Zellmembranen 

[0003] Da alle diese Eigenschaften aufgrund von genetischer Veranlagung, Gesundheitszustand oder Ein- 
fluss anderer Medikamente von Individuum zu Individuum schwanken konnen, und damit auch die Konzentra- 
tion des Wirkstoffs im Korper, kommt es auch bei der Wirkung von Medikamenten zu individuellen Unterschie- 
den. Diese konnen je nach Eigenschaft des Wirkstoffs unterschiedlich stark sein. Es ist bekannt, dass bei man- 
chen Wirkstoffen bei 50 % der Anwendungen keine therapeutische Wirkung erreicht werden kann. Die derzei- 
tige Vorgehensweise bei Beginn einer Therapie ist, die Verabreichung der Standarddosis und die darauffolgen- 
de Beobachtung des Patienten. Gegebenenfalls wird bei Nichteintritt eines Therapieerfolges versucht diesen 
durch schrittweise Erhohung der Dosis zu erreichen. Diese Vorgehensweise ist uneffektiv und kann zur Ge- 
fahrdung des Patienten fuhren. Letzteres gilt insbesondere auch fur Falle, in denen die Konzentration im Kor- 
per aufgrund individueller Faktoren wesentlich hoher ausfallt als im Normalfall, sodass es zu unerwunschten 
Nebenwirkungen kommt. 

[0004] Der Zusammenhang zwischen den individuellen Eigenschaften des Korpers und dem Verhalten eines 
pharmazeutischen Wirkstoff ist in vielen Fallen zumindest qualitativ bekannt. Besonderheiten mancher Ein- 
flussfaktoren, wie Korpergewicht oder Blutflussraten sind durch den behandelnden Arzt leicht diagnostizierbar 
(Gewicht) oder durch medizinisches Wissen, z.B. uber Anderungen der Durchblutung bei vorliegen einer 
Krankheit, abschatzbar. Dem Einfluss des Korpergewichts wird in manchen Fallen durch Verabreichung einer 
gewichtsspezifischen Dosis Rechnung getragen. Die derzeitigen Methoden lassen aber im Prinzip nur eine 
qualitative Berucksichtigung der individuellen Gegebenheiten zu. 

[0005] Im Fall des Einflusses von aktiven, biochemischen Prozessen ist die Lage wesentlich komplizierter. 
Hier kann es zu Veranderungen der Effektivitat der Prozesse durch unterschiedlich ausgepragtes Vorhan- 
densein der entsprechenden Proteine aufgrund genetischer Veranlagung, Krankheit oder externe Einflusse, 
wie z.B. andere Wirkstoffe kommen. Auch die Funktion der Proteine kann bei gleicher Expression z.B. durch 
genetisch bedingte Anderung der Proteinstruktur oder Wechselwirkung mit anderen Substanzen beeinflusst 
sein (vergl.: J. Licinio, M. Wong (Eds.), "Pharmacogenomics", Wiley-VCH, Weinheim 2002; A. D. Rodrigues, 
"Drug-Drug Interactions", Marcel Dekker, 2002). 

[0006] Die Wechselwirkung mit anderen Substanzen, insbesondere gleichzeitig verabreichten anderen Wirk- 
stoffen oder Nahrungsbestandteilen ist primar substanzabhangig und kann bei ausreichender Charakterisie- 
rung der Substanzen zumindest qualitativ vorhergesagt werden. Fur die wichtigsten Proteine existieren Listen 
von Substanzen, die deren Funktion entweder direkt oder durch Induktion oder Hemmung der Expression be- 
einflussen (vergl.: A. Schinkel, J. W. Jonker "Mammalian drug efflux transporters of the ATP binding cassette 
(ABC) family: an overview", Advanced Drug Delivery Reviews 55, 3-29 (2003); 1 Cytochrome P450 Drug In- 
teraction Table: http://medicine.iupui.edu/flockhart). In der Regel sind Einflusse auf den zu betrachtenden Wirk- 
stoff untersucht und in Form einer Kontraindikation in der Produktinformation vermerkt. 
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[0007] Schwieriger zu erfassen ist der Effekt durch die hinsichtlich der Struktur und Menge unterschiedliche 
Protein-Zusammensetzung der verschiedenen Organe eines bestimmten Patienten zum Zeitpunkt der Thera- 
pie. Hierzu musste im Prinzip die Expression eines bestimmten Proteins im betroffenen Organ bestimmt wer- 
den, was in der Regel nicht moglich ist, da keine entsprechenden Gewebeproben vorliegen. Es ist aber heute 
bekannt, dass Unterschiede in der Expression sowohl hinsichtlich der exprimierten Menge als auch der Prote- 
instruktur unter anderem auf Punktmutationen (Austausch einzelner Nukleotide in der genomischen DNA, so- 
genannte Single Nucleotide Polymorphisms „SNPs") zuruckzufuhren sind. SNPs konnen sich in den fur das 
Protein kodierenden oder in den die RNA-Transkription und damit auch die Expression des Proteins regulie- 
renden DNA-Sequenzen befinden. Daneben sind andere Arten von Mutationen (z.B. Insertionen, Deletionen) 
bekannt, die die Expression von Proteinen verandern konnen. Zudem ist bekannt, dass der genomischen DNA 
aufgepragte Methylierungsmusterdie Transkription/Expression verandern konnen. Derartige genomische Mar- 
ker wie z.B. SNPs lassen sich auf DNA-Ebene im Blut oder anderen Korperflussigkeiten nachweisen. In ver- 
schiedenen Fallen sind bereits speziell auf diesen Nachweis ausgelegte Verfahren, z.B. Biochips oder 
PCR-basierte Nachweisverfahren in der Entwicklung oder auf dem Markt. Dies eroffnet die Moglichkeit einer 
angepassten Dosierung auf Grundlage eines entsprechenden DNA-Tests direkt in der Arztpraxis oder als La- 
bormethode. Fur viele ADME-relevante Proteine liegen Informationen uberdie im Menschen auftretenden zu- 
gehorigen Gensequenzen und deren Modifikationen vor (vergl.: SNP-Datenbank im Internet verfugbar unter: 
http://www.ncbi.nlm.nih.gov/SNP/). Zu einem grofien Teil ist auch die Auswirkung der einzelnen Modifikationen 
auf die ADME-relevante Funktion bekannt (vergl.: R. G. Tirana, R. B. Kim "Pharmacogenomics of Drug Trans- 
porters" in J. Licinio, M. Wong (Eds.), "Pharmacogenomics", Wiley-VCH, Weinheim 2002). 

Aufgabenstellung 

[0008] Ein entscheidendes Problem bei der Bestimmung der optimalen individuellen Dosierung ist die gleich- 
zeitige, komplexe Abhangigkeit der intrakorporalen Konzentration von verschiedenen Einflussfaktoren. Wah- 
rend eine einzelne Abhangigkeit noch experimentell bestimmt und zum Beispiel in Tabellenform fur eine Do- 
sierungsentscheidung verwendet werden kann, ist dies bei mehreren, sich gegenseitig beeinflussenden Ab- 
hangigkeiten in der Regel hochstens qualitativ moglich. Diese Problematik lasst sich aber durch Verwendung 
einer computergestutzten Simulation zur Berechnung der Konzentrationen beheben. Eine hierfur geeignete 
Methode ist die sogenannte Physiologiebasierte Pharmakokinetische (PBPK) Simulation, mit der das Aufnah- 
me, Distributions-, Metabolisierungs- und Exkretionsverhalten (ADME) von Xenobiotika im Saugetierkorperauf 
Grundlage der physiologischen Voraussetzungen detailliert beschrieben werden kann. Fur einfache Fragestel- 
lungen, die nur passive Verteilungsprozesse im Korper berucksichtigen ist diese Methode schon lange bekannt 
und ausgiebig beschrieben (vergl.: G. E. Blakey, I. A. Nestorov, P. A. Arundel, L. A. Aarons, M. Rowland "Quan- 
titativ Structure-Pharmacokinetics Relationships: 1 . Development of a Whole Body Physiologically Based Mo- 
del to Characterize Changes in Pharmacokinetics Across a homologous Series of Barbiturates in the Rat", J. 
Pharmacokin. and Biopharm. 25, 277-312 (1997); R. Kawai, M. Lemaire, J.-L. Steiner, A. Bruelisauer, W. Nie- 
derberger, M. Rowland, "Physiologically Based Pharmacokinetic Study on a Cyclosporin Derivative, SDZ IMM 
125", J. Pharmacokin. and Biopharm. 22, 327-365 (1994)). Einige bekannte Modelle beschreiben zusatzlich 
den Metabolismus und auch den aktiven Transport in verschiedenen Organen als sattigbare, nicht dosislineare 
Prozesse. Ein Simulationsmodell, das durchgangig die Beschreibung solcher Prozesse in alien relevanten Or- 
ganen des Korpers berucksichtigtwird bei der PBPK-Simulationssoftware PK-Sim (Bayer Technology Services 
GmbH) verwendet. 

[0009] Die hier beschriebene Erfindung betrifft ein System aus der Kombination von einem Nachweissystem 
zur Bestimmung der ADME-relevanten genetischen Veranlagung des Patienten, einer PBPK/PD-Simulation 
und einer Datenbankfur Stoffeigenschaften (Fig. 1), das geeignet ist, dosisbezogen die individuelle Konzent- 
ration im Korper zu berechnen und aus dem Ergebnis die optimale individuelle Dosis vorzuschlagen. 

[0010] Gegenstand der Erfindung ist ein Verfahren zur Bestimmung der Dosierung von mindestens einem 
Wirkstoff auf Basis einer genetischen Analyse mit folgenden Schritten: 

a) Analyse (101) von spezifischen Gensequenzen, mittels eines gensequenzspezifischen Analysengerates, 
insbesondere eines sequenz-spezifischen Sensors, oder Bestimmung der Expression von Proteinen ent- 
weder uber RNA Transkription mittels quantitativer RNA-spezifischer Nachweismethoden oder direkte Mes- 
sung der Proteinexpression durch ein Protein-Analysegerat 

b) Zuordnung der Gensequenzen zu physiologischen Funktionen des menschlichen oder tierischen Kor- 
pers, insbesondere solchen physiologischen Funktionen, die Einfluss auf den Abbau, auf die Aufnahme, 
auf die Abgabe oder auf das Verteilungsverhalten des Wirkstoffes im Korper haben, 

c) Weitergabe der genetischen und der Zuordnungsdaten an einen physiologiebasiertes Pharmakokineti- 
sches Modell (PBPK-Modell) (108), 
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d) Eingabe von Wirkstoff-spezifischen Daten in das PBPK-Modell (108), 

e) Eingabe von charakteristischen Patientendaten, gegebenenfalls aus direkten Messungen am Korper 
oder in ein Rechnersystem (104), 

f) Berechnung von fur das PBPK-Modell notwendigen physiologischen Einflussparametern aus den Patien- 
tendaten unter Verwendung von Informationen aus physiologischen Wissensdatenbanken und Weitergabe 
der Parameter an das PBPK-Modell (108), 

g) Berechnung der individuellen Dosis aus den Daten nach den Schritten c), d) und f) unter Verwendung 
des PBPK-Modells (108). 

[0011] Durch die Gen-Test-Methode (101) wird die genetische Veranlagung des Patienten bezuglich der fur 
das ADME-Verhalten von Wirkstoffen wichtigen Gene oder Proteine bestimmt. 

[0012] Bevorzugt ist ein Verfahren, bei dem die Gensequenzen ausgewahlt sind aus solchen die die Proteine 
der Reihe: 

Metaboligsierende Enzyme, insbesondere Monooxygenasen der Cytochrome P 450 Familie, Phase II Enzyme, 
die polare Gruppen an die auszuscheidenden Molekuie anhangen, Aktive Transporter, insbesondere Multidrug 
Resistance Proteine z.B. P-Glycoprotein Familie oder Multidrug Resistance-Associated Proteins (MRP) oder 
Organic Anion Transporting Polypeptide Familie (OATP) oder Organic Anion Transporter Familie (OAT) oder 
Organic Cation Transporter Familie (OCT), oder Novel Organic Cation Transporter Familie (OCTN) oder Pep- 
tide Transporter Familie (PepT), oder Plasma Bindungs Proteine, insbesondere Serum Albumin und Glycopro- 
teine betreffen. 

[0013] Die wirkstoffspezifischen Daten sind besonders bevorzugt solche, die ausgewahlt sind aus der Reihe: 
Organ/Blut-Verteilungskoeffizienten, Membranpermeabilitat, kinetische Konstanten der Metabolisierungspro- 
zesse und/oder der aktiven Transportprozesse. 

[0014] Die charakteristischen Patientendaten gemaft Schritte) sind besonders bevorzugt ausgewahlt aus der 
Reihe: 

Korpergewicht, Korperoberflache, Korperfettgehalt, Alter oder Geschlecht, vom Normalzustand z.B. krank- 
heitsbedingt abweichende physiologische Funktionen wie beispielsweise renale oder hepatische Unterfunkti- 
on, Co-Medikation. 

[0015] Die physiologischen Einflussparameter gemaft Schrittf) sind besonders bevorzugt ausgewahlt aus der 
Reihe: 

Durchflussrate Q x von Blut durch das Organ X, Volumen V x des Organs X oder Permeabilitats-Oberfla- 
chen-Produkt (PxSA x ) fur das Organ X. 

[0016] Bevorzugt ist das PBPK-Modell ein Simulationsprogramm, das wenigstens die folgenden Funktionen 
simuliert: Darmaufnahme Bluttransport, Verteilung in Organe durch Permeation oder aktiven Transport, Meta- 
bolisierung, Ausscheidung uber Urin oder Galle. 

[0017] Gegenstand der Erfindung istfernereine Vorrichtung zur Bestimmung der Dosierung von Wirkstoffen, 
insbesondere nach dem erfindungsgemaUen oben beschriebenen Verfahren, aufweisend wenigstens ein gen- 
sequenz-spezifisches Analysengerat (101), eine hiermit verbundene Rechnereinheit mit einem Programm um- 
fassend ein pharmakokinetisches Modell (108), einer Wissensdatenbank (105) und Eingabemodulen (104) fur 
Patientendaten, dadurch gekennzeichnet, dass als pharmakokinetisches Modell (108) das PBPK-Modell(108) 
eingesetzt wird. 

[0018] Die wichtigsten ADME-relevanten Proteine sind: 

- Metabolisierende Enzyme: Monooxygenasen der Cytochrome P 450 Familie; Phase II Enzyme die polare 
Gruppen an die auszuscheidenden Molekuie anhangen. 

-Aktive Transporter: Multidrug Resistance (P-Glycoprotein Familie) (MDR), Multidrug Resistance-Associa- 
ted Proteins (MRP), Organic Anion Transporting Polypeptide Familie (OATP), Organic Anion Transporter 
Familie (OAT), Organic Cation Transporter Familie (OCT), Novel Organic Cation Transporter Familie 
(OCTN), Peptide Transporter Familie (PepT) 

- Plasma Bindungs Proteine: Insbesondere Serum Albumin und Glycoproteine 

[0019] Fur viele dieser Proteine sind Variationen in der genetischen Kodierung bekannt, die unterschiedlich 
haufig auftreten und unterschiedlich starken Einfluss auf die Funktion der Proteine und damit das ADME-Ver- 
halten eines Wirkstoffs haben. 
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[0020] Neben dem direkten Einfluss uber ADME-relevante Proteine istauch ein Einfluss durch genetisch be- 
dingte krankhafte Zustande moglich, die indirekt einen fur das ADME-Verhalten wichtigen Prozess beeinflus- 
sen. Auch solche genetischen Voraussetzungen konnen in die Gen-Untersuchung einbezogen werden, wenn 
der Zusammenhang zum Wirkstoffverhalten untersucht und im PBPK/PD-Modell beschreibbar ist. 

[0021] Die Gen-Test-Methode (101) selbst kann zum Beispiel eine Methode zurdirekten Bestimmung der Ex- 
pression der relevanten Proteine im Organgewebe, der Transkription relevanter RNA Molekule oder aber eine 
Methode zur Detektion von SNPs der DNAaus Proben von Korperflussigkeiten sein. Bevorzugt handelt essich 
dabei urn eine Biochip- oder PCR-basierte Anwendung. 

[0022] Die Ergebnisse des Gen-Test werden nach einer testspezifischen Methode ausgewertet (1 02) urn die 
benotigten Informationen uber die Beeinflussung von ADME-relevanten Prozessen zu erhalten. So wird ent- 
weder direkt das Expressions-Level der Proteine oder im Fall der DNA-Analyse uber bekannte Zusammenhan- 
ge die Auswirkung auf die Funktion oder Expression des entsprechenden Proteins bestimmt. Genomische 
Marker wie z.B. SNPs konnen auch zur Einteilung von Patienten in bestimmte Gruppen wie z.B. schnell oder 
langsam metabolisierende Patienten verwendet werden. Derzeit wird auch nach genomischen Markern ge- 
sucht, die die Einteilung von Patienten in Responder/lslon-Responder oder Patienten mit und ohne erwartete 
Nebenwirkungen bezogen auf bestimmte Medikamente oder Gruppen von Medikamenten ermoglichen. Der 
sich so ergebende Datensatz (103) wird dem PBPK/PD-Modell (108) als Eingangsdaten ubergeben. 

[0023] Weitere patientenspezifische Daten, die fur die Dosisberechnung relevant sind (104), sind manuell ein- 
zugeben. Bei diesen Daten handelt es sich urn Informationen, die sich durch Messung, Exploration oder Ana- 
mnese erheben lassen. Einige Bespiele sind: Korpergewicht oder Korperoberflache, Korperfettgehalt, Alter, 
Geschlecht. Die sich aufgrund dieser Daten ergebenden Parameterwerte des PBPK/PD-Modells werden in ei- 
nem nachfolgenden Schritt (106) unterZuhilfenahme einer Wissensdatenbank (105) uber die zugrunde liegen- 
den Zusammenhange berechnet. Diese Wissensdatenbank kann z.B. auch Informationen uber den Einfluss 
bestimmter Krankheiten auf ADME-relevante Prozesse enthalten. 

[0024] Eine mogliche Ausfuhrungsform des Moduls fur die manuelle Eingabe der Patientendaten konnte ein 
Eingabegerat mit menugefuhrte Nutzeroberflache sein, die dynamisch angepasst, abhanglg von den eingege- 
benen Information weitere benotigte Daten abfragt. 

[0025] Die fur die Simulation des ADME-Verhalten benotigten wirkstoff-spezifischen Daten fur das zu verab- 
reichende Medikament sind in einer weiteren Datenbank (107) abgelegt. Bei diesen Daten handelt es sich urn 
die im PBPK/PD-Modell enthaltenen, von den physiko- und biochemischen Eigenschaften des Wirkstoff ab- 
hangigen Parameterwerte. Diese wurden vorab direkt experimentell oder durch Anpassung des Simulations- 
modells an pharmakokinetische und/oder pharmakodynamische Daten bestimmt. Beispiele dieser Daten sind 
Organ/Blut-Verteilungskoeffizienten, Membranpermeabilitaten und die kinetischen Konstanten der Metaboli- 
sierungsprozesse und der aktiven Transportprozesse. 

[0026] Die zentrale Einheitdes Systems ist das PBPK/PD-Simulationsmodell mit dem die eigentliche Berech- 
nung der intrakorporalen Konzentrationen erfolgt. Der typische Aufbau eines PBPK-Modells ist in Fig. 2 ge- 
zeigt. Die prinzipielle Vorgehensweise ist, den Korper in einzelne Kompartimente zu unterteilen und den Aus- 
tausch von Wirkstoff-Substanz zwischen diesen Kompartimenten mit Hilfe von Massenerhaltungsgleichungen 
7) zu beschreiben. Als Kompartimente werden sinnvollerweise die einzelnen Organe gewahlt. Gegebenenfalls 
mussen auch Teile der Organe als Unterkompartimente definiert werden, wenn entweder der Stofftransport 
zwischen diesen limitiert sein kann, oder Informationen uber Konzentration separat gewonnen werden mus- 
sen. 

[0027] Bei diesen Massenerhaltungsgleichungen handelt es sich urn gewohnliche Differentialgleichungen. 
Typischerweise haben sie die Form: 

dC C 
V x —£- = Q,-C ar -Q x -£- (Gleichung 1) 



V x = Volumen des Organs X 

C x = Konzentration der Substanz in dem Organ X 

Q x = Durchflussrate von Blut durch das Organ X 

C ar = Konzentration der Substanz, die das Organ uber das arterielle Blut erreicht 
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K x = Verteilungskoeffizient der Substanz zwischen Blut und Organ X im Gleichgewichtszustand 

[0028] Sie beschreiben die Anderung der Konzentration im Organ X verursacht durch die mit dem Blutfluss 
Qx in das Organ transportierte Menge, die Verteilung zwischen Blut und Organgewebe, bestimmt durch den 
Verteilungskoeffizienten Kx und die mit dem Blutfluss wieder aus dem Organ abtransportierte Menge. 

[0029] Fur viele Pharmawirkstoffe ist die Verteilung in die einzelnen Organe dadurch limitiert, dass sie lang- 
samer durch die Zellmembranen permeieren als sie uber das Blut ins Organ transportiert werden. In diesem 
Fall sind die Organe in verschiedene durch Membranen voneinandergetrennte Unterkompartimente zu unter- 
teilen und man erhalt ein Modell entsprechend Fig. 3. Die zu berucksichtigenden Unterkompartimente sind 
Plasmavolumen (301), Rote Blutzellen (302), Interstitielles Volumen des Organgewebes (304) und Zellvolu- 
men des Organgewebes (306). Rote Blutzellen und Zellen des Organgewebes sind von Membranen (303), 
(305) umschlossen durch die die Wirkstoffmolekule permeieren miissen. Fur die durch Membranen umschlos- 
senen Kompartimente miissen in den Massenerhaltungsgleichungen fur den Stofftransport Permeationsterme 
nach dem 2. Fick'schen Gesetz berucksichtigt werden. Diese haben allgemein die Form: 

t^-tcc// f j-\ cell \ 

Vf ■ = PxSA x -\ C P J - (Gleichung 2) 



V x = Zellvolumen des Organs X 

CP" = Konzentration der Substanz im Blutplasma, die nicht gebunden ist 
C x e " = Konzentration der Substanz in den Zellen des Organs 

K x = Verteilungskoeffizient der Substanz zwischen Blut und Organ X im Gleichgewichtszustand 
PxSA x = Permeabilitats-Oberflachen-Produkt fur das Organ x 

[0030] Die aktiven Prozesse Metabolisierung und a ktiver Transport konnen zum Beispiel durch sogenannte 
Michaelis-Menten-Terme berucksichtigt werden, die die Kinetik der biochemischen Reaktionen beschreiben. 
Eine Einbeziehung des aktiven Transports setztein permeationslimitiertes Modell, wie oben beschrieben, vo- 
raus. Ein detailliertes Organmodell inklusive der aktiven Prozesse ist in Fig. 4 dargestellt. Einer oder mehrere 
Metabolisierungsprozesse (401) bewirken eine Abnahme der Konzentration der ursprunglichen Substanz. Die 
aktiven Transportprozesse (402), (403) werden so beschrieben, dass sie den Transport von Wirkstoffmoleku- 
len uber die Zellmembran, parallel zum passiven Permeationsprozess bewirken. Bei diesen Prozessen ist zu 
beachten, dass einwartsgerichtete (402) von auswartsgerichteten (403) zu unterscheiden sind. Entsprechend 
der Fig. 4 ist Gleichung 2 wie folgt zu modifizieren. 



(Gleichung 3) 

[0031] Dabeisind: 

— C pl = Michaelis-Menten-Term zur Beschreibung der Kinetik des 

(k m +Cf) x 

Einwartstransportes. 



Auswartstransportes . 



= Michaelis-Menten-Term zur Beschreibung der Kinetik des 
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— meto ™ x ■ C c ° = Michaelis-Menten-Tenn zur Beschreibung der Kinetik der 

(k + C x ) 

Metabolisierung 



V^ ax = Maximalrate des Prozesses y 

k y = Bindungskonstante des Wirkstoffs am Protein das den Prozess y bewirkt 
V x = Zellvolumen des Organs X 

CP" = Konzentration derSubstanz im Blutplasma, die nicht gebunden ist 
C x e " = Konzentration der Substanz in den Zellen des Organs 

K x = Verteilungskoeffizient der Substanz zwischen Blut und Organ X im Gleichgewichtszustand 
PxSA x = Permeabilitats-Oberflachen-Produkt fur das Organ x 

[0032] Nach dem gleichen Prinzip werden auch Organe mit spezielleren Funktionen beschrieben, wie z. B. 
der Magen-Darm-Trakt, die Niere, oder der Gallentrakt. Dabei sind in der Regel zusatzliche Parameter zu be- 
rucksichtigen, die die speziellen physiologischen Funktionen beschreiben. Im Fall des Darms ist zusatzlich die 
ortliche Variation von Grofien wie PxSA und pH-Wert des Darminhaltes zu berucksichtigen. 

[0033] Uberdie Konzentrationen Cx sind die Massenerhaltungsgleichungen des in Gleichungen 1-3 beispiel- 
haftdargestellten Typs fur die einzelnen Kompartimente und Unterkompartimente entsprechend dem Schema 
in Fig. 2 gegenseitig verknupft. So ergibt sich ein System von abhangigen Diffentialgleichungen in derZeit, das 
fur vorgegebene Anfangswerte numerisch gelost werden kann. Die Losungen dieses Gleichungssystems er- 
geben die Konzentrations-Zeit-Beziehungen fur alle im Modell enthaltenen Kompartimente. 

[0034] Zur Beschreibung einer pharmakologischen Wirkung kann weiterhin die Konzentrations-Zeit-Bezie- 
hung in dem Kompartiment, welches das biologische Target des Wirkstoffs enthalt, mit einem pharmakodyna- 
mischen Effekt verknupft werden. Typische Effektfunktionen sind z. B.: 

• Hyperbolische oder sigmoide Emax-Modelle: 

Effekt = E 0 + Emax C » 
EC*, + CI 

Effekt = pharmakologischer Wirkparameter (zeitabhangig) 

E 0 = Basiswert des pharmakologischen Wirkparameters 

E max = Maximalwert der pharmakologischen Wirkung 

EC 60 = Konzentration, bei der 50 % des maximalen Effekts erreicht sind 

C x = Konzentration am Wirkort (zeitabhangig) 

Y = Formparameter 

• Potenzfunktionen: Effekt = E 0 + (3 C x bzw. Log-Linear Modelle: 



Effekt = E 0 + p Ln(C x ) 



Effekt = pharmakologischer Wirkparameter (zeitabhangig) 

E 0 = Basiswert des pharmakologischen Wirkparameters 

P = Parameter fur die Steigung des Effekts als Funktion der Konzentration 

C x = Konzentration am Wirkort (zeitabhangig) 

Y = Formparameter 

• Wirkstoffinteraktionsmodelle wie z. B. partieller oder vollstandiger Antagonismus, etc. 

• Kombinationen der vorgenannten Modelle, mit denen z. B. multiple Wirkzentren, oder Rezeptor-Transdu- 
cer-Wechselwirkungen beschrieben werden konnen. 

[0035] Die Arbeitsweise des gesamten Systems zur individuellen Dosisberechnung ist nun folgende. Zu- 
nachst mussen die individuellen Werte der Parameter des PBPK/PD-Modells, die von der Physiologie oder 
Anatomie abhangen, bestimmt werden. Dazu werden die Ergebnisse des Gen-Tests (101) ausgewertet und 
die Proteine identifiziert bei denen bzgl. Expression oder Funktion Abweichungen gegenuber der Normalpopu- 
lation zu erwarten sind (102). Fur diese Proteine werden dann fur die betroffenen Organe uber bekannte und 
gespeicherte Relationen die Expression, bzw. die Effektivitat bestimmt und entsprechend die v max - und k m -Wer- 
te berechnet. Aus den sonstigen eingegebenen Daten des Patienten, z.B. Korpergewicht, Korperfettgehalt, 
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Geschlecht, Alter und ggf. Krankheitsbild werden unter Zuhilfenahme der zugehorigen in der Wissensdaten- 
bank (105) hinterlegten Relationen die Parameter V, Q, K, PxSA sowie weitere bei der Beschreibung spezieller 
Organe, wie Magen-Darm-Trakt, Niere etc. benotigte Parameter bestimmt. Dabei werden die in der Wirk- 
stoff-Datenbank (107) hinterlegten Standardwerte derwirkstoffspezifischen Parameter berucksichtigt, die dann 
entsprechend den individuellen Gegebenheiten moduliert werden. Gegebenenfalls ist auch der genetische 
bzw. krankheitsbedingte Einfluss auf Eigenschaften wie z.B. die Zusammensetzung der Zellmembranen und 
pH-Werte einzelner Kompartimente zu berucksichtigen, die ebenfalls Permeabilitaten PxSA und Verteilungs- 
koeffizienten K beeinflussen konnen. 

[0036] Ist der individuelle Parametersatz bestimmt, erfolgt die Simulation der Pharmakokinetik des zu verab- 
reichenden Wirkstoffs mit der Standarddosierung. Nach wirkstoff- und therapiezielabhangigen Regeln, die 
ebenfalls in der Wirkstoffdatenbank hinterlegt sind wird dann anhand der berechneten Konzentrationsverlaufe 
und ggf. des sich daraus ergebenden pharmakodynamischen Effekts entschieden, ob eine Anpassung der Do- 
sis notwendig ist. Ist dies der Fall wird eine geeignetere Dosis vorgeschlagen. Diese wird durch lineare Um- 
rechnung auf die im Korper zu erreichende optimale Dosis bestimmt, falls man sich im Regime der dosisline- 
aren Pharmakokinetik bzw. Pharmakodynamik befindet. Ist dies nicht der Fall wird durch automatische, schritt- 
weise Anderung der Dosis in der Simulation dies an das Optimum angepasst. Das Ergebnis dieser Optimie- 
rungen wird von System ausgegeben und kann dann zur Festlegung der Dosis verwendet werden. 

Erklarung der Figuren und Tabellen 

[0037] Fig. 1 : Prinzipieller Aufbau des Gesamten Systems zur Bestimmung der individuellen Dosierung. 

[0038] Fig. 2: Prinzipbild des Aufbaus des physiologiebasierten pharmakokinetischen (PBPK-) Simulations- 
modells Fig. 3: Zusammensetzung eines Organs im PBPK-Modell. 

[0039] Fig. 4: Prinzip der Beschreibung von aktiven Transportern und Metabolisierungsprozessen im PB- 
PK-Modell. 

[0040] Fig. 5: Simulierte Konzentrations-Zeitkurve (Linie) im Blutplasma von Patienten mit CC Polymorphis- 
mus in Exon 26 (C3435T) im Vergleich zu experimentell bestimmten Werten (Punkte). 

[0041] Fig. 6: Wie Fig. 5. Zusatzlich Konzentrationskurve fur Patienten mit TT Polymorphismus in Exon 26 
(C3435T) (graue Linie). 

Tabelle 1 : Experimented und durch Simulation berechnete Cmax-Werte fur Verabreichung 

von 0,25 mg Digoxin in Patienten mit CC und TT Polymorphismus in Exon 26 
(C3435T) 

Tabelle 2: Organvolumina und Blutflussraten 

Tabelle 3: Zusammensetzung der Organe entsprechend Fig. 3 

Tabelle 4: Substanzabhangige Parameter fur Digoxin 

Tabelle 5: Organspezifische substanzabhangige Parameter fur Digoxin 

Tabelle 6: Michaelis-Menten-Konstanten fur P-gp in der Darmwand 

Ausfuhrungsbeispiel 

[0042] Im Folgenden wird ein Beispiel behandelt, das zeigt, wie der Einfluss genetischer Veranlagung fur ak- 
tive Transporter auf die Pharmakokinetik von Wirkstoffen durch Simulationen beschrieben werden kann. 

[0043] Das am besten untersuchte und in der Literatur beschriebene Transportprotein ist p-Glycoprotein 
(P-Gp), das neben anderen Organen insbesondere auch im Darm exprimiert wird und dort Einfluss auf die Ab- 
sorption oral applizierter Wirkstoffe nehmen kann. Das zugehorige Gen wird in der Regel mit MDR1 benannt. 
MDR1 kann in verschiedenen Allelen vorliegen, wobei bekannt ist, dass diese zu unterschiedlicher Aktivitatdes 
zugehorigen Proteins fuhren (siehe Martin F. Fromm, The influence of MDR 1 polymorphisms on P-glycopro- 
tein expression and function in humans, Advanced Drug Delivery Reviews 54, 1295 (2002)). Fur P-Gp sind 
Tabellen von pharmazeutischen Wirkstoffen veroffentlicht, welche Substrat fur dieses sind (siehe z.B. C. J Ma- 
theny etal., Pharmacokinetic and Pharmacodynamic Implications of P-glycoprotein Modulation, Pharmacothe- 
rapy 21 , 778 (2001 )). Ein Beispiel fur ein solches Substrat ist Digoxin. Fur Digoxin ist bekannt, dass seine orale 
Absorption davon abhangt, welches Allel des MDR1 Gens vorliegt (siehe S. Hoffmeyer etal., Proc. Of the Na- 
tional Academy of Science USA, 97, 3473 (2000)). Konzentrations-Zeit-Kurven von Digoxin im Blutplasma fur 
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Patienten mit normal funktionaler MDR1-Sequenz sind in A. Johne et al., Clinical Pharmacology & Therapeu- 
tics, 66, 339 (1999) veroffentlicht. 

[0044] Es wurde ein PBPK-.Modell entsprechend Fig. 2 erstellt, mit dem. diese Konzentrationskurven durch 
Simulation beschrieben werden. In diesem Modell wird bei der Absorption der Substanz aus dem Darm neben 
der passiven Permeabilitat der Darmwand auch ein aktiver Transportprozess in Richtung des Darmlumens be- 
rucksichtigt. Unter Verwendung der in Tabellen 2-6 aufgelisteten Parameterwerte liefert das Simulationsmodell 
eine nahezu exakte Beschreibung der experimentell bestimmten Plasma-Konzentrationskurve (Fig. 5). 

[0045] Der in Tabellen 2-6 aufgefuhrte Parametersatz ist im Anwendungsfall in der Datenbank mit Wirkstoff- 
informationen (107) hinterlegt. Die Abschatzung des Einflusses von zum Beispiel abweichenden MDR1-Se- 
quenzen kann dann durch entsprechende Anderung derdadurch betroffenen Parameter des Simulationsmo- 
dells erreicht werden. Diese Anderung wird unter Berucksichtigung der Expressionsdaten (103) im Schritt „Pa- 
rameterbestimmung" (106) durchgefuhrt. 

[0046] Ein Beispiel fur einen Datensatz, in dem das Expressionslevel von P-Gp in der Darmwand in Abhan- 
gigkeit vom MDR1 Polymorphismus bestimmt wurde, ist in S. Hoffmeyer (2000) zu finden. Das Expressionsle- 
vel bestimmt wiederum die maximale Geschwindigkeit Vmax des Transportprozesses. Die Wissensdatenbank 
(105) enthalt also im Anwendungsfall die Zuordnung von relativen Vmax-Werten zu den Gensequenzen der 
verschiedenen Polymorphismen, die mit den substanzspezifischen absoluten Vmax-Werten aus der Daten- 
bank mit Wirkstoffinformationen (107) zu den im PBPK-Modell (108)zu verwendenden Parametern verrechnet 
werden. 

[0047] In S Hoffmeyer et al., (2000) sind zusatzlich Daten uber den Einfluss verschiedener Polymorphismen 
auf die Pharmakokinetik von Digoxin enthalten. So werden fur den Polymorphismus in Exon 26 (C3435T) die 
in Tabelle 1 aufgelisteten Cmax-Werte (Maximale Plasmakonzentration) angegeben. 

[0048] Ausgehend von der oben dargestellten Simulation des „Wildtypen" (C3435C) kann der Homozygote 
Typ TT durch Reduktion von Vmax urn 51 % entsprechend dem niedrigeren Expressionslevel simuliert werden. 
Das entsprechende Ergebnis wird in Fig. 6 dargestellt. Die sich in der Simulation ergebende Erhohung von 
Cmax urn 45 %, was der experimentell gefundenen Erhohung von 30 % gut entspricht (siehe Tabelle 1). 

[0049] Im Anwendungsfall wurden sich aus der Simulation des Typs TT ergebenden Cmax Werte mit in der 
Datenbank mit Wirkstoffinformationen (107) enthaltenen sicherheitskritischen Werten verglichen und gegebe- 
nenfalls eine reduzierte Dosierung vorgeschlagen. Urn den Dosierungsvorschlag festzulegen werden z.B. Si- 
mulationen mit iterativ veranderten Dosierungen durchgefuhrt, bis die pharmakokinetischen Kenngroften im si- 
cheren Bereich liegen. In Fallen, in denen sicher ist, dass eine lineare Abhangigkeit von der Dosis besteht kann 
der Dosierungsvorschlag auch durch linear Umrechnung bestimmt werden. 

Tabelle 1 



Polymorphis- 
mus 


Mittl. rel. Expres- 
sionslevel 


Cmax Experiment 
fug/11 


Cmax Simulaton 


CC (Wildtyp 


1275 


1,6 


1,72 


TT (Homozygot) 


627 


2,1 


2,49 
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Tabelle 2 





Voluroina 
fmll 


Blutflussratcn 
[ml/min] 


Venoser Blutpool 


250 


4670 


Arterieller Blutpool 


140 


4670 


Lunge 


670 


4670 


Magen 


150 


60 


Diinndarm 


640 


600 


Dickdarm 


370 


240 


Pankreas 


100 


60 


Milz 


180 


180 


Leber 


1710 


390 


Gallenblase 


20 




Niere 


720 


1133 


Gehirn 


1486 


700 


Herz 


330 


240 


Muskel 


30200 


550 




Organ 


Volumina 
[mil 


Blutflussraten 
[ml/min] 


Knochen 


12060 


167 


Haut 


3020 


50 


Fett 


10060 


300 


Hoden 


35 


2,6 
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Tabelle 3 



Volumenanteil 


f_vas 


fint 


fcell 


Fettgewebe 


0,010 


0,135 


0,855 


Gehirn 


0,037 


0,004 


0,959 


Magen-Darm-Trakt 


0,032 


0,100 


0,868 


Herz 


0,262 


0,100 


0,638 


Niere 


0,105 


0,200 


0,695 


Leber 


0,115 


0,163 


0,722 


Lunge 


0,626 


0,188 


0,186 


Muskel 


0,026 


0,120 


0,854 


Knochen 


0,041 


0,100 


0,859 


Haut 


0,019 


0,302 


0,679 


Bauchspeicheldriise 


0,180 


0,120 


0,700 


Milz 


0,282 


0,150 


0,568 


Hoden 


0,140 


0,069 


0,791 



Tabelle 4 



Parameter 


Wert 


Einheit 


Dosis 

Wasserldslichkeit 


0,25 
65 


mg 
mg/1 


inlrinsisische Leberclearance 


0,48 


ml/min/kg 


Freie Fraktion in Plasma 


0,73 




Passive Darmwand-Permeabi- 
litat 


1,4 10" 6 


cm/s 
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Tabelle 5 



Organ 


Organ/Plasma 
Verteilungs- 
koeffizient 


Permabilitats- 
Oberflachen- 
Produkt [ml/min] 


Magen 


8,34 


100 


Dunndarm 


8,34 


1000 


Dickdarm 


8,34 


500 


Pankreas 


10,57 


500 


Milz 


2,74 


250 


Leber 


9,35 


1000 


Niere 


7,19 


1000 


Lunge 


1,97 


0,3 


Gehirn 


14,21 


0,002 


Herz 


13,28 


1000 


Muskel 


2,33 


1000 


Knochen 


34,45 


28 


Haut 


13,49 


0,33 


Fett 


100,42 


14 


Hoden 


4,42 


3 








Tabelle 6 




Organ 


v 

Y max 
[ml/min] 


Km 
fmg/Il 




Duodenum 


0,086 


6,3 




Jejunum 


0,76 


6,3 




Ileum 


1,13 


6,3 





Patentanspriiche 



1 . Verfahren zur Bestimmung der Dosierung von mindestens einem Wirkstoff auf Basis einer genetischen 
Analyse mitfolgenden Schritten: 

a) Analyse (101) von spezifischen Gensequenzen, mittels eines gensequenz-spezifischen Analysengerates, 
insbesondere eines sequenzspezifischen Sensors, oder Bestimmung der Expression von Proteinen entweder 
uber RNA Transkription mittels quantitativer RNA-spezifischer Nachweismethoden oder direkte Messung der 
Proteinexpression durch ein Protein-Analysegerat 

b) Zuordnung der Gensequenzen zu physiologischen Funktionen des menschlichen oder tierischen Korpers, 
insbesondere solchen physiologischen Funktionen, die Einfluss auf den Abbau, auf die Aufnahme, auf die Ab- 
gabe oder auf das Verteilungsverhalten des Wirkstoffes im Korper haben, 

c) Weitergabe der genetischen und der Zuordnungsdaten an einen physiologie-basiertes Pharmakokineti- 
sches Modell (PBPK-Modell) (108), 

d) Eingabe von Wirkstoff-spezifischen Daten in das PBPK-Modell (108), 

e) Eingabe von charakteristischen Patientendaten, gegebenenfalls aus direkten Messungen am Korper, 
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f) Berechnung von fur das PBPK-Modell notwendigen physiologischen Einflussparametern aus den Patienten- 
daten unter Verwendung von in der Wissensdatenbank enthaltenen Informationen und Weitergabe der Para- 
meter an das PBPK-Modell (108), 

g) Berechnung der individuellen Dosis aus den Daten nach den Schritten c), d) und f) unter Verwendung des 
PBPK-Modells (108). 

2. Verfahren nach Anspruch 1, urch gekennzeichnet, dass die Gensequenzen ausgewahlt sind aus sol- 
chen die die Proteine der Reihe: 

Metabolisierende Enzyme, insbesondere Monooxygenasen der Cytochrome P 450 Familie, Phase II Enzyme, 
die polare Gruppen an die auszuscheidenden Molekule anhangen, Aktive Transporter, insbesondere Multidrug 
Resistance Proteine z.B. P-Glycoprotein Familie oder Multidrug Resistance-Associated Proteins (MRP) oder 
Organic Anion Transporting Polypeptide Familie (OATP) oder Organic Anion Transporter Familie (OAT) oder 
Organic Cation Transporter Familie (OCT) oder Novel Organic Cation Transporter Familie (OCTN) oder Pep- 
tide Transporer Familie (PepT), oder Plasma Bindungs Proteine, insbesondere Serum Albumin und Glycopro- 
teine 
betreffen. 

3. Verfahren nach einem der Anspruche 1 oder 2, dadurch gekennzeichnet, dass die wirkstoffspezifischen 
Daten solche sind, die ausgewahlt sind aus der Reihe: 

Organ/Blut-Verteilungskoeffizienten, Membranpermeabilitat, kinetische Konstanten der Metabolisierungspro- 
zesse und/oder der aktiven Transportprozesse. 

4. Verfahren nach einem der Anspruche 1 bis 3, dadurch gekennzeichnet, dass die charakteristischen Pa- 
tientendaten gemaR, Schritt e) ausgewahlt sind aus der Reihe: 

Korpergewicht, Korperoberflache, Korperfettgehalt, Alter oder Geschlecht. 

5. Verfahren nach einem der Anspruche 1 bis 4, dadurch gekennzeichnet, dass die physiologischen Ein- 
flussparameter gemaft Schritt f) ausgewahlt sind aus der Reihe: 

Durchflussrate Q x von Blut durch das Organ X, Volumen V x des Organs X oder Permeabilitats-Oberfla- 
chen-Produkt (PxSA x ) fur das Organ X. 

6. Verfahren nach einem der Anspruche 1 bis 5, dadurch gekennzeichnet, dass das PBPK-Modell ein Si- 
mulationsprogramm ist, daswenigstensdiefolgenden Funktionen simuliert: Darmaufnahme Bluttransport, Ver- 
teilung in Organe durch Permeation oder aktiven Transport, Metabolisierung, Ausscheidung uber Urin oder 
Galle. 

7. Vorrichtung zur Bestimmung der Dosierung von Wirkstoffen, insbesondere nach einem Verfahren ge- 
ma(3> einem der Anspruche 1 bis 6, aufweisend wenigstens ein gensequenz-spezifisches Analysengerat (101), 
eine hiermit verbundene Rechnereinheit mit einem Programm umfassend ein pharmakokinetisches Modell 
(108), einer Wissensdatenbank (105) und Eingabemodulen (104) fur Patientendaten, dadurch gekennzeich- 
net, dass als pharmakokinetisches Modell (108) das PBPK-Modell (108) eingesetzt wird. 

Es folgen 5 Blatt Zeichnungen 
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Anhangende Zeichnungen 
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Fig. 4 
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